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ABSTRACT 

We review five often used quad lens models, each of which has analytical solutions and 
can produce four images at most. Each lens model has two parameters, including one 
that describes the intensity of non-dimensional mass density, and the other one that 
describes the deviation from the circular lens. In our recent work, we have found that 
the cusp and the fold summations are not equal to 0, when a point source infinitely 
approaches a cusp or a fold from inner side of the caustic. Based on the magnification 
invariant theory, which states that the sum of signed magnifications of the total images 
of a given source is a constant, we calculate the cusp summations for the five lens 
models. We find that the cusp summations are always larger than 0 for source on the 
major cusps, while can be larger or smaller than 0 for source on the minor cusps. 
We also find that if these lenses tend to the circular lens, the major and minor cusp 
summations will have infinite values, and with positive and negative signs respectively. 
The cusp summations do not change significantly if the sources are slightly deviated 
from the cusps. In addition, through the magnification invariants, we also derive the 
analytical signed cusp relations on the axes for three lens models. We find that both on 
the major and the minor axes the larger the lenses deviated from the circular lens, the 
larger the signed cusp relations. The major cusp relations are usually larger than the 
absolute minor cusp relations, but for some lens models with very large deviation from 
circular lens, the minor cusp relations can be larger than the major cusp relations. 

Key words: gravitational lensing: strong - methods: analytical. 


1 INTRODUCTION 

Producing multiple images of distant quasars or galaxies by 
foreground galaxies or galaxy clusters is one of the most 
distinct qualities of strong gravitational lensing. For nonsin¬ 
gular le nses, it is w ell known that the total image number 
is odd jBur k<j 1 198ll l. If a point source lies within the cen¬ 
tral astroid caustic of the elliptic lens, there will be five 
images produced. There are two positive minima (minima 
point of Fermat potential, similar hereafter) images outside 
of the tangential critical curve, two negative saddle images 
inside of the tangential critical curve, and one positive max¬ 
ima image lying near the len s centre llBlandford fe Naravanl 
1 19861 : ISaha fe Willia ms 2003). However, the maxima image 
located near the lens centre is usually highly demagnified 
and faint, resulting in four observed images. 

There are some important magnification relations for 
the multiple image lenses. The magnification invariant 


means that, for some specific lens models, the sum of signed 
magnifications for all lensed i mages of a g iven point source 
is a constant, i.e., I = JT/A (lDala3ll998n . It is very inter¬ 
esting and surprising that the invariants are independent of 
most of the model parameters. For example, the magnifi¬ 
cation invariants of the point lens and Singular Isothermal 
Sphere (SIS) lens are 1 and 2 respectively, no matter how 
large the Einstein radii are and where the positions of the 
point sources are, as long as there are two images produced. 

The cusp and fold relations are local magnification 
relations compared with the magnification invariant. If a 
point source moves to the cusp from the inner side of 
the tangential caustic, three of the images will merge to¬ 
gether near the critical curve. The t hree close i mages have 
an as y mptotic magnification rela tion ( BlandforcUfc^lan^an 
19861: Schneider & Wcis^l 1 19921 : ISchneider. Ehlers fe Falco 


19921: iMad 119921 
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where p are the signed magnifications of the triple images 
A, B and C. Here, we define S C nsp and S'| cusp |, and name the 
numerator S cusp cusp summation, which will be frequently 
used in this work. If the point source infinitely approaches 
the cusp, the cusp relation R cusp will be close to 0. 

A similar magnification relation holds when the source 
lies near a fold caustic. In this case, two images lie closely 
together, straddling the critical curve. One of two images 
is a minima and the other one is a saddle. The fold 


(Blandford & Naravan! 

198f|: Schneider. Ehlers & Falco 

1992j; Mad 199^; 

Keeton. Gaudi & Petters 2005? 


Goldberg et al.[|20ld ) 


-Rfold = 


ftold 


Pa + 


S|foid| | A* a | + |me 


( 2 ) 


where fi are the signed magnifications of the double images 
A and B. Here, we define ftoid fold summation and another 
quantity S| fold | as before. If the source infinitely approaches 
the fold line, the fold relation i?f 0 id will also be close to 0. 

In some previous works, when the point source infinitely 
approaches the cusp or the fold, the numerators in E quations 
(1) and (2) are also considere d to be equal to 0 ( Zakharovl 
1 19951 : lAazami &; Pettersl 2009|^ Petters fe Wemerl 20ld 'l. In 


our recent work, we I Chu. Lin fe Yandl2013fl proved that 


Scusp and Sfoid are usually not equal to 0. Consequently, 
there are different signs in the numerators, so in the defi¬ 
nitions about the two relations we do not use the absolute 
value of the summed magnifications in the numerators as 
some other authors do. 

In strong gravitational lensing, the positions of most 
multiple images can be fitted adequately using simple 
smooth lens models. Neverthel ess, the observed flux ratios 
are more difficult to match (iKochanekI Il99ll b Actually, 
most of the observed fluxes of image pairs and triples 
disagree with the fold and cusp relations. The discrep¬ 
ancy between the predicted and observed flux ratios is 
commonly referred to as the anomalous flu x ratio problem 
llMao et al.l 200j^_ Congdon fe Keetonl 120051 : iMcKean et al.l 


l2007l : Shin fe Evansl 20081') . Currently the most favoured 


explanation of the flux ratio anomalies invokes the per- 


lensing galaxies (Mg 

10 & Schneider 199g|; Metcalf & Zhao 

2003; Kochanek & Dalai 

2004: Aazami & Natarajan 2006j: 

Maccio & Miranda 

2006 

; Chen. Kousliiappas & Zcntner 

20 111; Xu et all 200^, 

2015 

). In this work, based on the mag- 


nification invariant, we mainly study the cusp summation 
and cusp relation through five frequently used smooth quad 
lenses, and they may be helpful for our understanding of 
the anomalous flux ratio problem in another aspect. 


2 GENERAL REVIEW FOR THE FIVE QUAD 
LENSES 

We review five often used quad lenses in strong lensing, 
including Singular Isothermal Elliptical Density (SIED), 
Singular Isothermal Elliptical Potential (SIEP), Singular 
Isothermal Quadrupole (SIQ), SIS+shear, and Point-|-shear 
lenses. There are some similar properties for the five lenses. 
For each lens model, the radial critical curve degenerates 
into a point in the lens centre, and corresponds to the 
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<7 

Figure 1. The coefficients of multipole expansions of the SIED 
and SIEP lenses. For each panel, the seven curves describe the 
coefficients of modes m = 0, 2, 4, 6, 8, 10, 12 from top to bottom, 
respectively. 


pseudo-caustic dEvans fe Wilkinsonlll998l b Each of the five 
lenses has an astroid caustic which has four cusps and four 
folds, and each lens can produce four images at most for a 
single source. 

The detailed information about these lenses is shown 
in Table 1. Each lens model has two parameters, including 
one that describes the intensity of non-dimensional mass 
density, and the other one that describes the deviation from 
the circular lens. For the two SIE lenses, they are usually 
studied in the Cartesian coordinates, while for the last three 
lens models, it is more convenient to treat them in polar 
coordinates. 

The two SIE lenses are extended from the SIS lens. 
Here b is used as a constant parameter, which indicates 
the intensity of the mass density, q is the axial ratio of the 
SIE lens. The SIE lenses can be derived by changing 6 into 
\Jq 2 x 2 + y 2 through the SIS lens. For the SIED lens the 6 
was changed in the mass distribution k, while for the SIEP 
lens it was changed in the potential if. 

For the SIED and SIEP lenses, their convergence can 
also be written in the form of polar coordinates 


16 1 by/ 1 + e 

29 q 2 cos 2 4> + sin 2 (/> 2 6 y/1 — t cos 2<j> 
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Table 1. Five simple quad lenses. 


Lens model 

SIED“ (x,y) 

SIEP 6 {x,y) 

Parameters 

b, q (0 < q < 1) 

b, q (0 < q 7 1) 

Convergence 

r - b 

.. _ b q 2 (x 2 +y 2 ) 

2y/ q 2 x 2 +y 2 

2(q 2 x 2 +y 2 ) 2 / 2 

Deflection potential 

ip = xa x + ya y 

1 p = by/ q 2 x' 2 + y 2 

Deflection angle 

n - b tan -1 ( t 1-92 

bq 2 x 

aX yiZ^2 tan i 

^ y/q 2 X 2 -\-y 2 

Deflection angle 


n- - by 

° y ypv tann W q 1 ^+y ,lV ) 

" y/q 2 x 2 +y 2 

Shear 

7 = K 

7 = K 

Magnification 

/i _ 1 = 1 — 2 K, 

fi- 1 = 1 — 2k 

Critical curve 

q 2 x 2 + y 2 = b 2 

bq 2 (x 2 + y 2 ) _ , 

( q 2 x 2 +y 2 ) 3 / 2 

Major axis /3 CU s P 

b * tan-if^ 13 ^ 

1 y/l^/ 2 V 9 / 


Minor axis /3 CU sp 

-b + -J= tanh” 1 (y/T=q*) 

b — bq 2 

Naked cusp 

q < 0.39 

q < y/2/2 

Magnification invariant 

» 2.8 

2 


siq c (0,0) 

SlS+shear^ ( 6 , </>) 

Point+shear e ( 0 , 0 ) 

0 B , k (0 7 k 7 1) 

0 E, 7 (0 ^ 7 < 1 ) 

0 E, 7 (0 ^ 7 < 1 ) 

k= t|(1 + k cos 20) 

K 29 

K, = 7 r6^S(0) 

-0 = # E 0 — ^QftkOc os 2(f) 

if) = O^Q — cos 2 (f) 

i/j = 0 E In 6 — cos 2(f) 

“rad = 0 E ~ § 0 Ek COS 20 

Orad = 6 »e - 70 COS 20 

«rad = ~ cos 2< /> 

otan = §0 E fesin20 

atari = 70 sin 20 

atan = 70 sin 2 /> 

7 = K 

7 = ^ + 7 cos 2 (f) — i7 sin 2(f) 

7 = + 7 cos 2(f) — \ r y sin 2(f) 

a 

CN 

11 

1 

3 . 

/i —1 = 1 — 7 2 — ^(1+7 cos 2(f)) 

y ,- 1 = 1 - 7 2 - || - 2^7 cos 20 

6 = 0e H - #e k cos 2 (j) 

e 0 e 1+ 7 co 1 2 ^ 

(^■) 2 = \/l — 7 2 sin 2 2(f) — 7cos 2(f) 

/ 3 cusp = § 0 E fc 

/ 3 cusp = 20 E 

(flcusp - 0 E 7 Ty 7 

k > 0.6 

7> 1/3 

Never 

1 

2/(1 - 7 2 ) 

1/(1 -T 2 ) 


Q Kassiola & Kov 

ner (1993k Kormann. Schneider & Bartelmann (1994): Keeton & Kochanek (1998 

): Keeton. Mao &; Witt (2000|). 

^Kassiola & Kovner 

dl993h; IWitt, Mao &; Keeton (200olk c Kochanek j 19911); IWoldesenbet & Williams 

(2012ll: [Chu. Lin & Yand (i2013ll 

“Kovner ( 1987h : 

Finch et al. J2002I): Keeton. Gaudi &; Petterd (20031). TChans & Refsdal (19791. 19841): Schneider. Ehlers & Falco 


dl992h; [An fc Evan j ll2006h. -^Dalall lll99gi); IPaJal fc Rabinl d200lf) ; IWitt fe Mad d2000i') 


Einstein radius is determined by the monopole, the Ein¬ 
stein radius of the SIE lens is 0 E = no- For a given pa¬ 
rameter b, with increasing the ratio q, 0 E decreases for the 
SIED lens, while increases for the SIEP lens. Except the 
monopole do, with increasing q, all the coefficients of the 
two SIE lenses decrease. In addition, the coefficients of all 
the even modes decrease with increasing of t he mode m . 
The SIQ lens, also called SIS+elliptical lens (iDalall 1 19981 : 
IWoldesenbet fc Williamsll2012l l. can be thought as the low¬ 
est order multipole expansion of the two SIE lenses, and its 
coefficients only include the monopole ao and the quadrupole 
a 2 . 

In polar coordinates, the shear can be decomposed into 
two parts, tangent ia l or radi al shear 7 +, and skew shear -y x 
feernstein fc Nakaiimal 120091 1. We can write them into the 
complex form 7 = 7+ + ry x . The direction of the skew shear 
is rotated 45° from those of the tangential or radial shear. 
For the real part, when 7+ > 0, it is tangential shear, while 
oppositely it is radial shear. The radial shear is only obvious 
for the void or the lens with negative mass distributions. 
For example, the convergence k = cos 20/2(9, can produce 
tangential shear in the positive density region, and radial 


ksiep = 


bq 2 


1 by/1 + e( 1 — e) 


26 ( q 2 cos 2 0 + sin 2 0) 3 / 2 26 (1 — ecos 20) 3 / 2 

(4) 

Here, the parameter e is related to the axial ratio q by e = 
(l-g 2 )/(l + g 2 ). 

For each SIE lens, the convergence has the form k = 
Gsie( 0 )/ 20 . The shape function Gsie( 0 ) can be decom- 
posed i nto m ultipole s throug h Fourier transform method 
^Keeton, Gaudi fc Pettersll2003ll 


Gsie(0) = ^2 am cos mcf), 


m =0 


r 2n 


1 f 

a m = — Gsie( 0 ) cos rrujxicf). 
^ 7r Jo 


(5) 


Based on the symmetry of the SIE lenses, the coefficients 
a-m of the odd modes are equal to 0. Therefore, there are 
only even modes in the expansions, and their phases 0 m are 
all equal to 0 . 

Figure 1 shows the Fourier expansion coefficients of 
the shape functions of the SIE lenses. The strength of the 
monopole is the largest one among all modes. Since the 
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shear in the negative region. In addition, the external shear 
7 in the last two lens models can bring both 7 + = 7 cos 2 <j> 
and 7 x = —7 sin 2 </ components. 

For any lens with a convergence of the form of G(4>)/29, 
one can decompose it into multipoles, similar to Equation 
(5) . For each mo de of this lens, except m = 1, it has 7 = k 
IlChu et al.l l2013h . These shears are all tangential or radial 
ones 7 +, not including skew shear 7 X . According to the su¬ 
perposition principle, we can conclude that, as long as the 
function G((/>) do not include the m = 1 mode, the shear of 
the lens is 7 = k, and the magnificat ion is/r = 1/(1 — 2k). 
In fact, IWitt. Mao fc Keetonl (|200Cll l found that the lens 
potential with if> = 9F(rf>) form has the magnification of 
fx = 1/(1 — 2k). It is consistent with our conclusion, because 
for a convergence G{(j>)/29 with no monopol e, the lens po¬ 
tential can be written in the form ijj = 9F(tf>) llEvans fe Wittl 
l200ll. l2003h . 

The last two lenses are derived by adding a uniform ex¬ 
ternal shear on the SIS or the point lens. The poin t +sliear 
lens i s also called Chang-Refsdal lens dChang fe Refsdall 
Il979l . ll984lh There is usually a minus sign before 7 in few 
other studies. We use it in this form to let the major axis 
lie along the A-axis, and minor axis along the T-axis. For a 
general strong lens, when the 9 is infinite large, the magnifi¬ 
cation n should have a positive sign. Therefore, it needs the 
external shear 7 to be smaller than 1. However, it cou l d also 
b e st udied with 7 > 1 for extreme cases, as in lAn fe Evand 

(l2006h . 

In addition, the convergence ti9^5{6) of the point lens 
or the point+shear lens is derived through the relation of 
the two-dimensional Dirac delta function V 2 In \6\ = 2x 18(d). 
Here for the point lens, we do not consider the Schwarzschild 
radius of the poin t mass, near which the deflection angles 
are more complex (IVirbhadra fe Ellis 2 000 ll . 

From the deflection angles and critical curves one can 
easily calculate the angular distance /3 C us P from the cusp 
to the source centre, and the angular distance /3 pse u from 
pseudo-caustic to the source centre on the axes. When the 
two quantities are equal to each other, we can derive the 
critical value of the second parameter for naked cusp ap¬ 
pearing. 

The lens equation (3 = 6 — a includes two independent 
equations. It describes the transformation between the lens 
plane (9,cj)) and the source plane (/ 3,ip ), and can also be 
written in polar coordinates as 

/3 2 = (9 - Qrad ) 2 + a?an, ( 6 ) 

tan(</> - ip) = Qtan (7) 

v C^rad 

dChu et al.l l2013h . Using these lens equations in the Carte¬ 
sian coordinates or the polar coordinates, together with the 
critical curves, one can also calculate the caustics for the 
SIE lenses or the other three lenses. 


3 THE CUSP SUMMATIONS FOR SOURCE 
ON THE CUSPS 

As we know, when a point source is exactly on the cusp, the 
cusp summation is usually not equal to 0. However, it is diffi¬ 
cult to derive the magnifications of the triple images related 


to cusp summation using traditional methods. Fortunately, 
we have the magnification invariants for the five lens mod¬ 
els. Therefore, we can derive the cusp summation through 
the differences between magnification invariant and the fi¬ 
nite magnification of the fourth image. The magnification 
invariants of the SIEP, SI Q, SIS+shea r, Poi nt+shear lenses 
have been calculated by iDalall (1 1998) and iDalal fe Rabinl 
d200lj) . The magnif ication invariant of SIED lens is given by 
IWitt fe Mad (120001) . These magnification invariants are only 
valid when four images are produced. 

For the SIED lens, the magnification invariant is 
llWitt fe Madl2000l ~) 


Isied = 


1 — u /tan 


1 — v /tanh 1 v 


• 2 . 8 , ( 8 ) 


where u = -y/2e/(l — e), and v = \j2tj (1 + e). The magnifi¬ 
cation invariant slightly depends on the q, and also slightly 
depends on the position of the source. 

For a point source on the cusp, it has one image with 
finite magnification, and three images merged together with 
infinite magnifications. After we derive the finite magnifi¬ 
cation of the point source, together with the magnification 
invariant, we can get the cusp summations of the triplets 
both on the major and minor cusps 


2 + u/ tan 1 u 2 

2 — 2u/ tan -1 u 1 — v/ tanh -1 v ’ 
2 2 + v/ tanh -1 v 

1 —u/tan -1 w 2 — 2 v/tanh -1 v 


( 9 ) 


When e ~ 0.73 or q « 0.39, the cusp summation of the major 
cusp is equal to the magnification invariant. When e ss 0.37 
or q « 0 . 68 , the cusp summation of the minor cusp is equal 
to 0 . 

Similarly, for the SIEP lens, we can also derive the cusp 
summations of the triplets on the major and minor cusps 


3-2 q 2 
2 — 2q 2 
2-3 q 2 
2-2 g 2 


I (5 + 7 } ’ 

iO- 


( 10 ) 


When e = 0.2 or q = %/6/3, the cusp summation of the minor 
cusp is equal to 0. For the SIQ lens, the cusp summation of 
the triplet on the major/minor cusp is 

Scuap = g(l ± ^)- ( 11 ) 

For the SIS+shear lens, the cusp summation of the triplet 
on the major/minor cusp is 


S 


cusp - 


57 + 1 
47(1 - 7 2 )’ 


( 12 ) 


When 7 = 0.2, the cusp summation of the minor cusp is 
equal to 0. For the Point+shear lens, the cusp summation 
of the triplet on the major/minor cusp is 


= 4 7 ±(l-7 2 ) 

87(1 + 7 ) 


(13) 


When 7 ~ 0.24, the cusp summation of the minor cusp is 
equal to 0. When 7 approaches 1, the minor cusp summation 
will be close to 0.25. 

For each of these lenses, the cusp summation does not 
depend on the first parameter of the lens, and only depends 
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Figure 2. The red and green curves show the cusp summations on the major and minor cusps, respectively. The black curves mean the 
magnification invariants. 


on the second parameter which describes the deviation from 
the circular lens. Figure 2 shows the cusp summations on 
the cusps for the five lens models based on Equations (9)- 
(13). Apparently, the cusp summations S m s P are usually not 
equal to 0 , and are small quantities compared to the infinite 
magnifications of the triple images. Therefore, the value of 
Scusp is very easily to be ignored as higher order infinitesi- 
mals, when treat with the magnifications the three images 
llSchneider fe Weisslll992l ; I Zakh arovll 199, il l . 

If these lenses tend to the circular lens, the major and 
minor cusp summations will have infinite values. The cusp 
summations of the major cusps are always larger than 0 . 
When it is smaller than the magnification invariant, it means 
naked cusps appear. For the Point+shear lens, naked cusp 
will never appear, because the pseudo-caustic of this lens 
model is a circle with infinite large radius. In addition, unlike 
the major ones, the cusp summations of the minor cusps can 
be larger or smaller than 0 except for the SIQ lens. 


4 THE CUSP SUMMATIONS AND 

RELATIONS FOR SOURCE ON THE AXES 

Through the difference between magnification invariant and 
the magnification of the fourth image, one can also calculate 
the cusp summation on the major or the minor axis. For each 
lens model, we define 2 = /3//3 C us P both on the major and 
minor axes. Because the analytical solutions for the SIED 
and point+shear lenses are very complex, we do not study 
them here. 


For the SIEP lens, the cusp summations on the major 
and minor axes are 


Scusp — 1 + 


(1 - q 2 ){l + z)’ 
2 


5cusp_1 (1 — g 2 )(l + s)' 


(14) 


For the SIQ lens, the cusp summations on the major and 
minor axes are 


3(fc + 1) 
4fc(l + z) ’ 
3 (k - 1 ) 
4fc(l + z) 


(15) 


For the SIS+shear lens, the cusp summations on the major 
and minor axes are 


_ 7(3 + 2z) + 1 

27(1 - 7 2 )(1 + «)’ 
7(3 + 2z) - 1 
27(1 — 7 2 )(1 + z)' 


(16) 


Figure 3 shows the cusp summations of the three lens 
models based on Equations (14)-(16). We can find that the 
cusp summations on the cusps are the smallest ones for 
source on the major axes, while the cusp summations on the 
cusps are the largest ones for source on the minor axes. The 
cusp summations do not change very much if the sources are 
slightly deviated from the cusps. 

If a point source lies on the axes, one can easily calcu¬ 
late the positions and the magnifications of the two images 
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Figure 3. The solid and dashed curves show the cusp summations 
on the major and minor axes, respectively. 


locating on the axes of the lens plane. However, through 
traditional methods, it is difficult to calculate the magnifica¬ 
tions of the other two images which have same magnification 
values. Now, since we know the magnification invariants for 
these lens models, we can also analytically derive the mag¬ 
nifications of the two images lying off the axes. 

After deriving the magnifications of the two side images, 
we can also calculate the cusp relations on the major and 
minor axes. For the SIEP lens, the cusp relations on the 


major and minor axes are 


Rcmsp 2 + q 2 + z + z 2 — q 2 z 2 


Rcusp - 


V 


— 1 — 2 q 2 — q 2 z + z 2 — q 2 z 2 


- 1 , 

+ 1 . 


( 17 ) 


For the SIQ lens, the cusp relations on the major and minor 
axes are 


3(1 + fc)(l — z) 

9 + k + 3z + 3 kz + 8 kz 2 ’ 
3(1 — fc)(l — z) 

—9 + k — 3z + 3kz + 8 kz 2 


(18) 


For the SIS+shear lens, the cusp relations on the major and 
minor axes are 


R 


cusp 


R 


cusp 


4 + 47 

3 + 7 + Z + 72-I- 27 z 2 
4 — 47 

3 + 7 — z + 72 + 27a 2 


- 1 , 

+ 1 . 


(19) 


Figure 4 shows the cusp relations of the three lens mod¬ 
els based on Equations (17)-(19). For each of the three lens 
models, the cusp relation also does not depend on the first 
parameter of the lens, and only depends on the second pa¬ 
rameter. For each lens, major cusp relations approach 0 only 
from positive values. However, for the minor cusp relation, 
they can approach 0 from positive or negative v alues. In fact, 
the early work lKeeton, Gaudi fe Pettersl J2003I') also implied 
that the situation for minor cusp is not certain, and regret¬ 
fully they studied the absolute value of the R CU s P . Neverthe¬ 
less, the minor cusp relations of the SIQ lens only approaches 
0 from negative values. 

As shown in Figures 3 and 4, on the major axes, the 
changes of S CU s P and R cusp against the second parameters 
are usually opposite. E.g., in Figure 3(a), the solid red curve 
is the lowest one, while in Figure 4(a), the solid red curve 
is the highest one. On the minor axes, the changes of S CU s P 
and -Rcusp are similar. It has been proved that, when e, k, 
and 7 are close to 0, the cusp summations S cusp are infinitely 
large. Through Equations (17)-(19), we find even these three 
parameters equal to 0, as long as z = 1, the cusp relations 
Rcusp are still equal to 0. It means that when these lenses 
tend to circular lens, 5| cusp | are higher order infinities com¬ 
pared to the Scusp. 

In observations, for the major cusp relation with a posi¬ 
tive sign, the flux summation of the two side images is larger 
than that of the middle image. For the minor cusp relation 
with a positive sign, the flux of the middle image is larger 
than the sum value of two side images, while with a nega¬ 
tive sign, the sum flux of the two side images is larger than 
that of the middle image, which is similar to the major cusp 
in observation. From Figure 4, we find that the larger the 
lens deviated from circular lens, the larger the signed cusp 
relations, for both on the major and minor axes. The minor 
cusp relations are more sensitive to the second parameters 
than the major cusp relations, especially for the SIEP and 
SIS+shear lenses. The minor cusp relations of these two lens 
models can be larger than 0, and can be even larger than 
the major cusp relations. 

The discriminant for triple images in being a major or a 
minor cusp type can be found in the distances of the images 
from the lens centre. If the distance from the triplet to lens 
centre is larger than that from the singlet, it is a major cusp 
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Figure 4. The solid and dashed curves show the cusp relations 
on the major and minor axes, respectively. 


type, otherwise it is a minor cusp type d Met calf fc Amaral 
1 20121) . RXJ0911+0551 feade et al.ll 19971: iBurud et al.l llS 


is the only cusp type lens system whose source lying near 
the minor cusp, and the signed cusp r ela tion of this lens 
syste m is R cuap = —0.192 ±0.011 llKeeton. Gaudi fe Pettersl 
l200l the negative sign is given by us). If it can be fitted 
by the SIEP lens, we can conclude that the axes ratio q is 
very large. There is no sample in which the middle image 
has larger flux than the sum of two side images observed by 
now. If it can be find, it can be confirmed as minor cusp 
type with a larger deviation from circular lens. 


5 CONCLUSIONS AND DISCUSSION 

The four-image lens systems are very important and 
are very common in the observations of lensed quasars 
llRusin fc Tegmarkll200ll ; Idaeskens fe Surdeill2002fl . We re¬ 
view five quad lens models, each of which has analytical 
solutions and can produce four images at most. Each of the 
five lenses has two parameters, including one that describes 
the intensity of the mass distribution, and the other one 
that describes the deviation from the circular lens. Using 
the magnification invariants of these lens models, we calcu¬ 
late the cusp summations for the five lenses. We find that 
for a point source on the cusp, the cusp summation is al¬ 
ways larger than 0 for the major cusp, while can be larger 
or smaller than 0 for the minor cusp. If these lenses tend 
to the circular lens, the major and minor cusp summations 
will have infinite values, and with positive and negative signs 
respectively. 

In this study, we calculate the cusp summations on the 
axes for SIEP, SIQ, SIS±shear lenses, and find that the cusp 
summations on the cusps are the smallest ones for sources on 
the major axes, while the cusp summations on the cusps are 
the largest ones for sources on the minor axes. The cusp sum¬ 
mations do not change very much if the sources are slightly 
deviated from the cusps. 

In addition, through the magnification invariants, we 
also calculate the magnifications of the two side images for 
source on the axes, and then derive the signed cusp rela¬ 
tions on the axes. We find both on the major and the minor 
axes that the more the lenses deviated from the circular 
lens, the larger the signed cusp relations. When the point 
source moves to the minor cusp infinitely, the minor cusp 
relation has two ways to approach 0, i.e., from positive or 
from negative value. The changes of Sens P and R cusp against 
the second parameter are usually opposite for sources on the 
major axes. When these lenses tend to circular lens, 5| cusp 
are higher order infinities compared to the 5 CU sp. 

The analytical results show that, the major cusp rela¬ 
tions are usually larger than the minor cusp relations, but 
for some lens models with larger deviation from circular lens, 
the minor cusp relation can be larger than the major cusp 
relation. In some previous numerical work, the major cusp 
relations are much more easily larger than the minor cusp 


relations (Bradac et al. 

2004; Amara et al. 2000; Xu et al.1 

20091; Metcalf & Amara 

20121). We guess that in most of the 


numerical simulations the cllipticity of projected haloes is 
too small to let the minor cusp relation to be much larger. 


The cusp summation and relation do not depend on 
the first parameter of the lens. Therefore, if we change the 
redshifts of the lens body or the source (which is equivalent 
to multiply a constant on the first parameter), it usually 
does not change the cusp summation or the cusp relation. 
We guess the fold summation and fold relation are similar 
to those properties of the cusp. It should be noted that, 
all the conclusions in this work are based on these simple 
lens models, and we do not consider the nonsingular core or 
substructures in the real lens bodies. There is no doubt that 
they can influence the cusp or fold magnification relations 
significantly. 

In future, based on the magnification invariants, we can 
also calculate the fold summations of these quad lenses for 
source exactly on the fold lines, through the finite magnifi- 
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cations of the other two images which can be derived using 
numerical method. We expect the summed magnification of 
fold image pair to change continuously along the fold line. 
Furthermore, closer to the major cusp, the fold summations 
will be smaller, while closer to the minor cusp, the sum¬ 
mations will be larger, and the fold summations will have 
infinite values with negative and positive signs near the m a- 
jor and minor cusps respectively (iChu, Lin fe Yanal2013h . 

Based on the cusp summations of the five lenses, we can 
bring forward a question: is the major cusp summation al¬ 
ways larger than 0 for any type of major cusps? Or in other 
word, is it independent of the lens models? This is a very 
interesting mathematical problem. Solving this question can 
help us to understand the singularity theory or caustic meta¬ 
morphoses in strong gravitational lensing. In addition, it 
is also very interesting to study the magnification summa¬ 
tions for a poi nt source on the higher order singul a rities of 
the caustics llPetters. Levine fe Wambs ganss 200ll; IWernerl 
120091 : lAazami. Petters fc Rabinll201ll l. In these cases, more 
than three images with infinite magnifications will merge 
together, and the dependence of the magnification summa¬ 
tions on the lens parameters are more complex than those 
of the cusp. 
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